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Abstract 


The  model  we  describe^has  the  same  general  features  of  the  PILOT 

dynamic  macro-economic  model  of  U.S.  designed  to  assess  the  long  term 

impact  of  foreign  competition,  innovation,  modernization,  and  energy 

needs.  We  derive  the  aggregate  demand  function  of  final  consumer  from 

individual  demand  functions  in  order  to  state  its  mathematical  properties; 

we  then  estimate  its  parameters  by  a  fit  to  empirical  data.  The 

equilibrium  conditions  are  those  of  the  Arrow-Debreu  model,  the  only 

unusual  feature  is  that  investors  calculate  their  rate  of  return  using 

discounted  normalized  prices  of  future  periods.  If  investors  choose  to 

normalize  intra-period  prices  in  the  usual  way  by  requiring  that  they  sum 

to  unity  (or  equivalently  their  average  value  is  unity),  the  inverse  demand 

functions  turn  out  to  be  non-integrable .  Equally  satisfactory  from  the 

investors'  point  of  view,  is  for  them  to  choose  instead  to  normalize 

_  _  _  i/2 

intra-period  prices, by  making  these  equal  to  n/(n'Hn)  where  H 

is  a  given  positive-definite  matrix  and  u  is  the  vector  of  intra-period 

'  *  11  Au  v. 

prices.  In  the  latter  case,  we— show  that  the  inverse-demand  functions  are 
integrable  and  derive  a  utility  function  for  the  economy  which  if  maximized 
subject  to  the  physical-flow  constraints  implies  the  equilibrium 


conditions. 
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We  begin  by  describing  the  structure  of  the  t  =  1,  . ..,  T  time 
period  model.  It  has  the  same  general  feature  as  the  PILOT  model  of  U.S. 
See  references  [12,  13]. 


EQUILIBRIUM  MODEL.  For  t  =  1,  ...»  T: 
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where  variables  with  a  hat  over  them  are  resp.  the  slack  vectors  for  the 
inequalities  (l),  ...,  (4)  and  5  t+^  is  the  discount  factor.  Relations 
(5)  are  the  complementary-slackness  conditions.  See  references  [1,4,14]. 
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Based  on  joint  work  with  P.H. 


McAllister  and  J.C. 


Stone. 


The  matrix  structure  is  skew  symmetric  except  for  the  term  in  lower 
right  diagonal.  It  is  of  interest  to  note  that  if  the  economy  were  driven 
by  a  utility  function  of  the  form  U  =  E  6  Ut(Xt>,  t*ie  Kuhn-Tucker 
conditions  derived  by  maximizing  U  subject  to  primal  physical-flow 
conditions  (1)  and  (2)  would  have  exactly  this  structure  where  5Ut/5Xt 
*  F^CX^).  If  the  latter  condtions  hold,  we  say  the  vector  functions 
Ft(Xt)  in  the  model  are  integrable.  See  references  [3,10]. 

(1)  states  that  the  capacity  requirements  to  meet  production  and 
investment  levels  Y  in  period  t  must  not  exceed  capacity  carried  down 
from  period  t-1  plus  any  exogeneous  capacity  supplied  k^.  For  period  1, 
the  term  +DqYq  *8  omitted;  k^  is  the  initial  capacity  (endowment) 
vector. 

(2)  states  that  the  consumption  vector  Xf  of  the  final  consumers 
must  not  exceed  the  net-output  of  production  after  investment.  The 
consumption  vector  of  government  services  is  treated  as  part  of  AfcYt  and 
is  not  shown  separately. 

(3)  states,  in  case  of  an  intra-period  production  activity  j,  that 

production  level  Y t ( j )  will  rise  to  the  point  of  non-profitability  and  if 

strictly  non-prof itable  will  not  be  used  because  of  the  complementary 

slackness  condition  Y  Y  *  0.  In  the  case  of  an  inter-period  investment 

activity,  the  investor  must  receive  his  discounted  rate  of  return  5  or  he 

won't  invest.  Prices  it  are  normalized  intra-period  t  prices 
-t+1  T 

discounted  by  6  .On  the  right  hand  side,  is  the  discounted 


value  of  future  capacity  (endowments).  For  t  -  T  this  term  is  set  -  0. 
Note:  Undiscounted  prices  will  be  denoted  by  itt* 

(4)  is  an  equation  for  Xfc  >  0.  It  relates  prices  to  consumption. 

The  direct  demand  function  Is  an  expression  of  the  form  Xt  **  ^(i^,  n£Xt) 
where  Sf  Is  a  homogeneous  function  In  of  degree  0  and  n£Xt  is 

the  attained  level  of  aggregate  income  1^.  It  follows  from  the  homogenlty 
of  Sf  that  the  inverse  demand  function  F(X)  (expressing  it  as  a  func¬ 

tion  of  X^)  can  only  determine  prices  within  a  scale  factor.  It  is  the 
freedom  to  select  the  scale  factor  which  allows  us  to  choose  the  formula 
for  normalizing  intra-period  prices  «t.  We  define  F(X)  as  equal  to 
normalized  intra-period  prices. 

We  assume  that  the  Investors  will  want  to  use  discounted  normalized 
intra-period  t  prices  in  calculating  their  rate  of  return.  If 
intra-period  prices  are  normalized  by  scaling  them  so  that  their 

average  is  unity  then  the  average  of  the  discounted  prices  it  is  5  t+*. 
Since  it^  -  6  t+*  Ft(Xt),  when  Xt  >  0,  this  implies  that  F^  satisfies 
eFt(Xt>  =  I  for  all  choices  of  Xfc  where  e  *  (1/n,  ...,  1/n)  and  n  is 
the  number  of  components  of  Xt. 

If  this  usual  way  of  normalizing  intra-period  prices  is  used,  i.e.,  so 
that  their  average  is  1/n  or  as  we  prefer  1,  it  will  turn  out,  however, 
that  the  equilibrium  problem  (1),  ...,  (5)  is  non-integrable,  i.e.,  there 
does  not  exist  a  utility  function  which  if  maximized  subject  to  physical- 
flow  conditions  (1)  and  (2)  and  (X^Y^)  0  yields  the  equilibrium  solu¬ 

tion.  We  will  show,  on  the  other  hand,  that  there  is  another  way  to  do 
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the  normalization,  equally  satisfactory  from  the  view  point  of  the 
Investor,  that  is  integrable  and  from  which  it  is  easy  to  derive  an 
aggregate  utility  function  for  the  economy. 

Derivation  of  the  aggregate  demand  function 

Let  lt  be  the  value  of  endowments  used  to  produce  the  consumption 

vector  in  period  t.  In  terms  of  the  prices  of  the  model,  the  value 

T  T 

of  endowments  available  to  period  t  is  o  (D  ,Y  .  +  k  )  *  a  B  Y 

w  L  X  t  [  C 

T 

and  the  value  passed  down  to  period  t+1  is  at+jD  Y  .  Their  difference 

is  by  definition.  By  (1),  (2),  and  (5),  it  is  clear  1^  -  ntXt* 

aggregate  take-home  income.  We  will  often  omit  the  subscript  t  to 
simplify  the  notation;  thus  we  write  I  -  1^. 

We  will  use  the  index  i  to  denote  the  1-th  consumer;  I*  *  I*  to 

denote  his  personal  income  in  period  t;  and  X1  -  X*  to  denote  his 

consumption  vector.  Let  U^(X^)  denote  his  utility  function,  and  let 
be  his  assumed  constant  share  of  aggregate  income  I.  The  price  vector 
corresponding  to  the  components  of  X*  is  denoted  by  n  *  x  . 

The  budget  constraints  are 

(6)  n'X1  -  o^I  ; 

T.  w’X1  -  I  ,  where  *  1. 

At  equilibrium  we  have  for  individual  i, 

(7)  5 (u1(x1)  -  x^’x1)]  /  ax1  -  0  , 

where  the  Lagrange  multiplier  X  is  chosen  so  that  his  budget  constraint 
x'X1  »  Ojl  is  satisfied. 
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We  approximate  U*(X*)  in  the  neighborhood  of  an  equilibrium  solution 
by  a  general  quadratic  function  of  the  form 

(8)  U1(Xi)  -  (M1S1)T  X1  -  (1/2)  (xV  M1(X1)  +  Constanti  , 

where  M*  in  general  is  a  symmetric  non-singular  matrix.  We  assume  M* 
to  be  positive  definite  and  S*  to  be  a  fixed  vector  such  that  X*  <  <  S* 
for  any  X*  attainable  by  the  model.  Of  course  S*  *  S*  is  strictly 
positive;  it  could  be  different  for  different  t.  In  the  PILOT  model 
-  ji  S*  grows  proportional  to  population  size,  i.e.,  as  the  number  of 
individuals  i  grows . 

We  substitute  the  approximation  (8)  into  (7)  obtaining 

(9)  M1(S1-X1)  •  X»x  , 

(10)  S^X1  i  X * H^x  ,  where  H1  -  (M1)”1  . 

Note  that  H*  »  (M*)  *  is  also  symmetric  and  positive  definite. 

We  can  now  use  the  budget  constraint  to  determine  X.  Multiplying 

(10)  by  it'  on  the  left,  setting  x'X1  »  a^I,  we  can  solve  for  X  and 
substitute  into  (9).  This  yields  the  local  approximation  for  the  demand 
function  of  Individual  i  in  the  neighborhood  of  the  equilibrium  solution 
as  a  function  of  prices  and  aggregate  Income  I: 

.  .  w'S1  -  a. I  . 

(11)  S1  -  X1  -  ( - — — )  H1*  . 

x  'H  x 
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Summing  over  all  1  and  setting  S  *  £S*  and  X  =  EX*  yields  the 
local  approximation  to  the  aggregate  demand  function. 

.  it’CS1  -  a  S)  a 

(12)  S  -  X  =■  [J  - j  “  n  ]x  +  (x'S-I)  •  — j  H  jit  . 

i  x'H  x  i  x'H  x 

Denoting  the  bracket  terms  above  by  G  and  H,  we  have 

(13)  S  -  X  *  Git  +  (it '  S  -  I)  •  Hit  . 

Note  G,  H  are  square  symmetric  matrices  because  they  are  weighted  sums 
of  square  symmetric  matrices  H*  whose  elements  are  functions  of  it 
only.  It  is  easy  to  verify  that  G,  H  have  the  following  properties: 

(14)  it ' Git  =  0  ,  it  'Hit  =  1 

for  all  x.  Moreover,  H  is  positive  definite  because  It  is  a  non-negative 
sum  of  positive-def inlte  symmetric  matrices  H*.  The  individual  demand 
functions  (11)  and  the  aggregate  demand  function  (13),  for  fixed  prices  it, 
are  locally  linear  in  I*  -  at^I  and  I  resp. 

We  now  make  a  fit  to  empirical  data  to  see  if  the  local  linear 
approximations  can  be  extended  to  a  broad  range  of  I.  For  this  purpose 
we  will  need  survey  data  of  personal  consumption  as  a  function  of  take-home 
Income  at  fixed  prices  and  we  will  need  the  distribution  of  take-home 
income  a^.  See  references  [2,  7,  8]. 
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Personal  consumption 

At  fixed  survey  year  prices  x^  *  (x^,...  x^****  *  persona-*- 

consumption  C^(I^)  during  the  survey  year  of  item  k,  say  food  measured 
In  physical  units,  is  known  from  survey  data  as  a  function  of  personal 
income  I*,  see  Figure  1: 


FIGURE  1:  Consumption  of  Food  as  a  Function_of  Personal  Take-home  Income 

(prices  are  fixed  at  x  -  x^) ) 

Since  the  sum  of  the  consumption  of  item  k  weighted  by  prices  x^  2  ® 
over  all  items  adds  up  to  take-home  income,  it  follows  that  if  some  curves, 
like  food,  display  a  decreasing  slope  with  increasing  income  I*  then 
others  must  show  an  increasing  slope  with  increasing  income.  In  PILOT, 
expenditure  patterns  0^(1*)  of  individuals  for  some  future  period  t  at 
constant  prices  Xq  are  assumed  to  be  the  same  as  the  survey  year  with 
some  adjustments  to  reflect  any  known  trends  in  "taste”. 
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Take-home  income 


The  distribution  of  take-home  income  I*  expressed  in  base-year 
dollars  has  changed  over  the  years  in  U.S.  but  generally  retains  the  same 
shape  except  it  spreads  out  proportionally  as  per-capita  take-home  income 
1  increases.  That  is  to  say,  the  proportional  share  of  endowments 
has  remained  more  or  less  constant  over  the  years,  [8,2]. 


FIGURE  2:  Distribution  of  Personal  Income  I*  =  Ij  and  I 2 
when  average  income  is  1^  and  resp. 

Letting  p( I | I )  denote  the  distribution  of  personal  take-home 
income  when  per  capita  income  is  I,  we  are  assuming  that  p(ljl)  can 
be  derived  from  p< I | Iq )  by 

(15)  p( I | X)  =  0- PC9 I l IQ) 

where  9  -  Iq/I. 


9 


(16)  X  (I)  -  /  p( X J X)  C  (I)  dl 

1=0 

00 

=  e  /  p<ex| x0)  c^i)  dl  ,  e  =  X0/I  . 

1=0 

Note  that  X  (I)  is  per-capita  consumption  as  a  function  of  per-capita 
income  I  expressed  in  survey-year  dollars  and  assumes  prices  are  fixed 
at  Xq.  Period  t  prices,  expressed  in  survey  year  dollars,  may  differ 
from  Xg.  Later  on  we  derive  how  it  varies  with 

At  fixed  survey-year  prices  Xg,  using  personal  consumption  data  as 
a  function  of  take-home  income  and  the  observed  distribution  of  income,  the 
functions  Xt^(l)  have  been  computed  in  the  manner  described  for  over 
a  hundred  commodities  by  M.  Avriel  and  for  aggregated  key  commodities 
(using  more  recent  data)  by  P.H.  McAllister.  These  functions  turn  out  to 
be  remarkably  linear.  See  Figure  3  and  references  [7,2,8]. 
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I  *  per  capita  take-home  Income 


FIGURE  3:  Consumption  of  Food  as  a  Function  of  Per  Capita 
Take-home  Income  (prices  are  fixed  at  it  *  x^). 

Global  versus  Local  Fit 

At  fixed  survey  year  prices  itg,  individual  demand  functions  do  not 
appear  to  be  linear  with  I1,  see  Figure  1.  However,  we  obtain  an 
excellent  global  fit  to  the  aggregate  demand  functions  using  a  linear 
function  of  per  capita  income  I,  see  Figure  3;  or  multiplying  1  by 
population  size  to  obtain  I,  by  a  linear  function  in  I.  Therefore  we 
accept  (13)  as  the  form  of  our  global  fit  to  the  aggregate  demand 
function.  We  think  of  S  -  St  as  a  kind  of  "satiation"  vector  many  times 
larger  than  any  X  *■  Xfc  attained  in  any  period. 

Returning  now  to  our  aggregate  demand  function  (13),  we  now  postulate 
that  as  aggregate  income  I  approaches  it ' S ,  the  income  sufficient  to 
purchase  the  satiation  vector  S,  aggregate  consumption  X  tends  to  S. 


! 

In  (13),  if  we  fix  prices  x  and  let  I  ♦  x'S,  we  observe  that  Git  =  0 
i  for  all  x  so  that  we  can  drop  the  first  term  of  (13).  Now  H,  in  the 

second  term,  is  a  symmetric  positive-definite  matrix  whose  elements  depend 
on  x  with  the  property  that  x'Hx  =  1  for  all  x;  it  follows  that  the 
!  general  form  of  H  is 

(18)  5  ■  » 

I 

where  H  is  a  positive-definite  matrix  whose  elements  can  depend  on  x. 

2 

Our  global  fit  to  the  aggregate  demand  function  thus  reduces  to 

(19)  S-X  «  •  Hx  ,  X  <  <  S  . 

H  HTC 

For  purposes  of  estimation  of  S  and  H,  we  assume  H  to  be  a  constant 
1  matrix  as  well  as  symmetric  and  positive  definite. 

To  obtain  the  inverse  demand  function  F(X)  that  expresses  x  in 
terms  of  X,  we  solve  (19)  for  x.  Since  I  -  x'X,  the  right-hand  side  is 
a  homogeneous  function  in  x  of  degree  0  implying  that  x  can  only  be 
determined  within  a  scale  factor.  Clearly  x  is  proportional  M(S-X) 
where  M  *  H  *  and  we  are  free  to  choose  the  proportionality  factor  so 
that  the  prices  of  the  model  are  automatically  normalized  before  discount¬ 
ing.  In  the  model  relation  (4),  x  »  6  t+1  F(X)  when  X  >  0  so  that  F(X) 
denotes  normalized  intra-period  prices. 


2 

McAllister  has  estimated  S  and  H  in  (19)  using  22  years  of 
empirical  data  of  per-capita  consumption  patterns,  per-capita  take-home 
income,  and  prices  under  the  assumption  that  H  is  a  constant  matrix  that 
is  symmetric  and  positive  definite.  See  reference  [3]. 


If  investors  calculate  their  rate  of  return  based  on  intra-period 
prices  itt  normalized  by  it fc/  e«t  where  e  -  (1/n,  1/n),  then 

the  scale  factor  is  chosen  so  that  normalized  intra-period  prices  F(X) 
satisfy 

(20)  F(X)  -  eM(-s_'x)  *  M(S-X)  ,  M  -  H_1  . 

Note  eF(X)  =  1  for  all  X  and  F(X)  does  not  depend  on  the  scaling  of 
H.  Note  that  when  prices  it*  *  (I,  1,  . 1)  that  (1,  1,  . ..,  1) 
are  normalized  prices. 

However,  if  Investors  calculate  their  rate  of  return  based  on  intra- 

-  -  -  1/2 

period  prices  n t  normalized  by  flc/(n^  H  *t)  ,  the  scale 

factor  is  chosen  so  that  normalized  intra-period  prices  F(X)  satisfy 

(21)  F(X)  -  - ^ - rTT  •  M(S-X)  ,  M  -  H_1  . 

[(s-xr  M(S-X)]1/Z 

Note  that  [F(X)  H  F(X)]  '  =  1  for  all  X  but  that  F(X)  does  depend 
on  scaling  of  M.  We  can  rescale  H  and  hence  M  in  (21)  so  that  if 
intra-period  prices  it*  ■  (1,  1),  they  will  also  be  (1,  1) 

after  normalization,  i.e.,  satisfy  In*  H  n*J  “  1.  Therefore  H  is 
rescaled  so  that 

(22)  Rescaled  H  -  H/(n*  H  it*) 
and  M  -  [Rescaled  H]  *  is  then  used  in  (21). 
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Both  ways  to  normalize  prices  appear  to  be  equally  satisfactory  from 


the  view  point  of  the  investor  figuring  his  rate  of  return.  However  if 
(20)  is  used,  we  will  give  an  easy  proof  below  that  no  utility  function  for 
the  economy  exists;  whereas  the  Interesting  thing  is  that  if  (21)  is  used, 
there  is  one. 


Proof:  Assume,  on  the  contrary,  that  a  utility  function  does  exist  for 

(20).  Consider  a  one  period  model  so  that  we  are  maximizing  the  utility 
U(X)  subject  to  the  primal  constraints  (1)  and  (2).  Further  suppose 
has  only  two  components  so  that  X^  -  (X^,  Xi2^*  I'et  S1  *  ^11’  ^12^’ 

S1-X1  *  (^li'Xij,  ^12~ ^12^*  11 1  m  ^n*  ^  m  fmi j  1  is  a  2*2 

symmetric  non-singular  matrix.  Let 


(23) 


X  - 


V1  *  Su  “  X11  »  v2  "  S12  "  X12  * 

At  a  maximum  the  Kuhn-Tucker  conditions  5U/3X 
S-V,  we  have  from  (20) 


V  -  (Vj,  v2)  . 

■  x  hold.  Since 


(24)  au/axu  -  xn  -  (muV1  +  m12V2)/D  , 

(25)  au/axi2  -  wn  -  (m12V1  +  m22V2)/D  , 

where  the  denominator  D  *  ^mn+mi2^Vl  +  ^ml2+m22^V2* 

In  order  for  a  utility  function  to  exist,  the  second  partial 

2  2 
a  U/axnax12  computed  from  (24)  should  agree  with  3  U/ax^dX^  computed 

from  (25)  for  all  choices  of  Xj^,  X^.  Setting  these  2nd  partials  equal 


to  each  other,  we  obtain 


12  mllVl  +  mj2V2)  (ffl12  *  m22) 
D  D2 


m12  (ml2Vl  +  m22V2)(mil  +  ml2) 
D  D2 


2  2 

which  reduces  to  (“nm22  ”  “12^  ^Vi  “  V2^°  E  °*  whlch  does  not  hold 
for  all  choices  of  (X^,  X^)  because  (vl_v2)  “  (S^-Xj  1^“^S12~X12^ ’  a 
contradiction. 


On  the  other  hand,  if  we  normalize  intra-period  prices  x  by 

-  -  -  1/2 

x/(x'Hx)  ,  then  the  equilibrium  problem  (1),  (5)  is  equivalent  to 

solving  the  convex-programming  problem: 


■MATHEMATICAL  PROGRAMMING  MODEL. 


Find  Minimum  -U(X) ,  (Xt,  Y  )  >  0: 


(27)  -U(X)  -  l  6"t+1[(S  -X  )T  M  (S  -  X  )] 

t-1  t  c  r  r  c 

subject  to  primal  flow  conditions  for  t  *  1,  ...,  T: 


Corresp. 


+  Vt 
-  Vt 


i  +  Dt-lYt-l  +  kt 

<  0 


ot  >  ° 

n  >0 
t  — 


Because  M  are  positive  definite  matrices,  it  is  not  difficult  to 
prove  each  term  of  (27)  is  a  convex  function  in  and  hence  their  sum 
-U  is  also.  References  [5,  6,  9,  11 J  discuss  the  existence  of  solutions 


for  convex  programs  and  the  techniques  for  their  solution.  The  primal  and 


dual  variables  of  the  optimum  solution  satisfy  the  Kuhn-Tucker  conditions 
which  are  precisely  those  of  the  equilibrium  problem  (1),  . (5). 

We  conclude  that  the  economy  will  grow  If  It  has  the  technology  and 
Initial  endowments  to  grow  and  If  It  pays  to  trade  off  movement  of  the 
consumption  vector  Xt  away  from  the  "satiation’’  vector  St  in  earlier 
periods  for  considerably  larger  movements  towards  the  satiation  In  later 
periods  where  the  measure  of  disutility  function  for  period  t  is  given  by 

(29)  "W  "  S’t+1l<St-Xt)T  Mt(St-Xt)]1/2  . 

By  (4)  and  (21),  -Ut(Xt)  -  1tt^t_Xt^  80  C^8t  t*'e  disutility  *8  the 
discounted  additional  aggregate  income  needed  to  purchase  the  "satiation" 
vector  —  i.e.,  the  more  additional  income  required  the  lower  the  ’’standard 
of  living”. 
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The  model  we  describe  has  the  same  general  features  of  the  PILOT 

dynamic  macro— economic  model  of  U.S.  designed  to  assess  the  long  term  impact 

of  foreign  competition,  innovation,  modernization,  and  energy  needs.  We 

derive  the  aggregate  demand  function  of  final  consumer  from  individual 

demand  functions  in  order  to  state  its  mathematical  properties;  we  then 

estimate  its  parameters  by  a  fit  to  empirical  data.  The  equilibrium 

conditions  are  those  of  the  Arrow-Debreu  model,  the  only  unusual  feature  is 

that  investors  calculate  their  rate  of  return  using  discounted  normalized 

prices  of  future  periods.  If  investors  choose  to  normalize  intra-period 

prices  in  the  usual  way  by  requiring  that  they  sum  to  unity  (or  equivalently 

their  average  value  is  unity),  the  inverse  demand  functions  turn  out  to  be 

non-lntegrable.  Equally  satisfactory  from  the  investors'  point  of  view,  is 

for  them  to  choose  instead  to  normalize  intra-period  prices  by  making  these 
-  -  -  1/2 

equal  to  it/(it'Hn)  where  H  is  a  given  positive-definite  matrix  and 
x  is  the  vector  of  intra-period  prices.  In  the  latter  case,  we  show  that 
the  inverse-demand  functions  are  integrable  and  derive  a  utility  function 
for  the  economy  which  if  maximized  subject  to  the  physical-flow  constraints 
implies  the  equilibrium  conditions. 
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